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Abstract 

The problem of the Kosterlitz-Thouless (KT) transition in the highly frus- 
trated XY Kagome antiferromagnet is solved. The problem is mapped onto 
that of the KT transition in the XY ferromagnet on the hexagonal lattice. The 
transition temperature is found. It is shown that the spin correlation function 
exponentially decays with distance even in the low-temperature phase, in con- 
trast to the order parameter correlation function, which decays algebraically 
with distance. 
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Generally, XY spins on two-dimensional lattices undergo a Kosterlitz-Thouless (KT) 
transition As a rule, physics of this transition does not depend on details of lattice 

structure. In the low temperature phase, pairs of spin direction singularities, vortices, with 
opposite topological charges form quasi- molecules. Above the transition temperature, the 
quasi-molecules decay into a vortex gas. In the low-temperature phase, the spin correlations 
decay algebraically with distance, above the transition temperature, they decay exponen- 
tially. Formally, the universality of the KT transition follows from the possibility to describe 
low-energy states of two-dimensional XY spin systems in terms of the nonlinear cr-model. 
The XY antiferromagnet on the two-dimensional Kagome lattice is an exception from this 
class. It has infinitely many ground states, therefore its low-temperature properties cannot 
be described by the nonlinear a-model. This makes the problem of a KT-like transition 
in the XY Kagome antiferromagnet a special one, which is substantially more complicated 
than the standard theory of the KT transition. The very possibility of a KT transition in 
such an unusual system is under question. 

The problem of the KT transition in the Kagome antiferromagnet was first addressed 
by Huse and Rutenberg ||. They suggested that the order parameter for the KT transi- 
tion is e 3te where 9 is the angle of a spin. This order parameter is invariant with respect 
to any arbitrary choice of ground states, which are a subset of local 27r/3 spin rotations. 
Therefore this order parameter can change smoothly in the plane even though spins rotate 
locally at multiple 27r/3 angles. An indirect evidence of the KT transition in the Kagome 
antiferromagnet was obtained from MC simulation J|. 

A network of Josephson junctions with the 7r-phase shift can be mapped onto the anti- 
ferromagnetic XY model as well. Experimental studies of artificial networks of Josephson 
junctions on the Kagome lattice make this problem especially appealing ||. 

In this paper, we examine the KT-like transition in the XY Kagome antiferromagnet. 
We find that the KT transition in the Kagome antiferromagnet does exist, we evaluate the 
transition temperature, and we show that the spin correlation function behaves itself in an 
unusual way. 
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FIGURES 




FIG. 1. The Kagome lattice (filled dots) with antiferromagnetic bonds (continuous lines) and 
the dual lattice (circles) and its bonds (dashed lines). 

The common treatment of the KT transition in the continuous limit |6||^] does not account 
for the structure and high degeneracy of the system and therefore it is not applicable in this 
case. In order to take into account the special structure of the Kagome lattice we follow the 
lattice approach developed by Jose, Kadanoff, Kirkpatrick, and Nelson (JKKN) M. The 
Kagome lattice consists of triangles and hexagons (fig. 1). The Hamiltonian of the Kagome 
antiferromagnet can be represented as a sum of squares of the total spins in triangles of the 
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nearest neighbors, Sa: 

#=1/2 J£(Sa) 2 , (1) 

A 

where A numbers the triangles of the nearest neighbors (fig. 1). Each spin participates in 
two triangles. The ground state energy is equal to zero and there are infinitely many ground 
states with Sa = 0. In any ground state, the angles between neighboring spins are equal to 
±27r/3. Vortices are structural defects in the ground state spin pattern with kernels localized 
in the centers of hexagons or triangles. 

Our plan is as follows. Following JKKN, we represent the partition function in terms of 
integer-valued currents on the lattice and introduce the dual lattice to the Kagome lattice, 
with sites located in the centers of triangles and hexagons. A very important step is to 
exclude vortices located in the centers of triangles from further consideration because the 
energy of such vortices is approximately Jin (2 a/3) higher than that of vortices residing on 
hexagons. The neglect of such vortices makes sense because, as we show further (Eq. ([12])), 
the KT transition temperature, T c = 0.0756J is small compared to the energy of a vortex 
residing on a triangle, and the contribution of such vortices is of order (2\/3) ~ 7 • 
10~ 8 <C 1. This small parameter allows one to integrate out currents in triangles in the 
partition function. After this step, only currents in hexagons and new variables, chiralities, 
located on triangles and equal to ±1 are left. The chirality distinguishes the clockwise and 
counter-clockwise configurations of spins in each triangle. We show that summation over 
chiralities renormalizes the partition function expressed in terms of currents in hexagons, 
however it does not change the low-temperature partition function drastically. We find that 
the degeneracy of ground states results in a decrease of the transition temperature compared 
to that of the XY ferromagnet on the hexagonal lattice. The order parameter correlations 
below the KT transition temperature decay as a power of distance. We show that in the 
Kagome antiferromagnet the spin-spin correlations decay exponentially, in contrast with 
those in "usual", non-degenerate XY magnets. The reason for such a fast decay of the spin- 
spin correlations is that the symmetry of the pair spin correlation function is lower than 
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that of the order parameter and those correlations are sensitive to the multiplicity of ground 
states, that are mainly disordered. 

The partition function of the XY Kagome antiferromagnet can be represented as a sum 
over the lattice's bonds 

Z((3) = f e-^^ cos[eir) - 9{r+a)] Y[de(r), (2) 

J r 

where r marks positions on the Kagome lattice, a are the three lattice vectors directed 
along the antiferromagnetic bonds between nearest neighbors, 6 r are the spin angles, and 
(3 = JS 2 /2T is the dimensionless inverse temperature. The 2ir periodicity of the angle 
variables allows one to expand Z(/3) in Fourier series with the coefficients I n (r,a)( — fl), where 
I n (x) is the modified Bessel function and integer numbers n(r, a) are located on bonds 
connecting nearest neighbors r and r + a. 

One can integrate over the angles 9(r) in the partition function. This results in the 
following representation for the partition function: 

m= E n4(r,a)(-/W E n(r,a)), (3) 

{n(r,a)} (r,a) \a=r— r' j 

where n(r — a, —a) = —n(r, a). The numbers n(r, a) obey the conservation condition at each 
site of the lattice: 

5>(r,a) = 0. (4) 

a 

Following JKKN || , we introduce the dual lattice to the Kagome lattice in order to account 
for the conservation conditions ([|). The sites of the dual lattice, R, are located at crossings 
of lines perpendicular to the bonds that connect nearest neighbors on the Kagome lattice 
and cross the bonds in their middle points (fig. 1). This construction provides a one-to-one 
correspondence between bonds of the Kagome lattice and bonds of the dual lattice. The 
conservation conditions (f|) can be resolved if one introduces integer-valued currents J(R) 
circulating in each triangle and in each hexagon of the the Kagome lattice. The currents can 
be assigned to sites R of the dual lattice located in the centers of triangles and hexagons. 
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A current along a bond (r, a), n(r, a), is equal to the sum of currents in one triangle and 
in one hexagon that share the bond (r, a). The currents J(R) obey the Kirchhoff's rule: 
the sum of currents arriving at each site is equal to the sum of departing ones. Thus, the 
conservation conditions (^) are fulfilled through the Kirchhoff's rules. 

The summation over currents {n(r, a)} with the conservation laws (|j) in the partition 
function ([D is equivalent to summation over currents {J(R)} located on sites of the dual 
lattice 

Z (P) = E II im+A )+ j(K)(-P), (5) 

{J(R)} (R,A) 

where A are the lattice vectors of the dual lattice. 

From this point, our way deviates from that by JKKN. The product II(r,a) over the 
bonds can be factorized as a product over all triangles of the Kagome lattice and triple 
products of the Bessel's functions corresponding to the bonds in each triangle because each 
bond of the Kagome lattice belongs to one and only one triangle. This allows one to represent 
the partition function as follows: 

z(p)= e n e n ^(R t+ A h , )+ w-/?)- (6) 

{J(R h )} (R t ) J(R t ) h'=l 

Here we separate the sums over hexagon and triangle currents, J(Rh) and J(Rt), with 
centers Rh and Rt, and h' numbers three hexagons surrounding each triangle Rt. 

We consider sums of the triple products of the Bessel functions that appear in Eq. (||) 
and recall that I n (—/3) is n-th Fourier harmonic of e _/3cos<? \ In the low-temperature limit 
(P ^> 1), the integrand in each triple product has extrema at 4>h = 2ncr h /3, a h = ±1 
(h = 1,2,3). Thus, one arrives at the following asymptotic formula 
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e n / j(R t +A h ,)+j(R t )(-/?) ~ 

J(Rt) h'=l 

_ / ^TT^Rt) * 

2^ exp i ^ ^R t +A h , 

<r(R t )=±l I ° h'=l 

— tjJ^-^j^O^)). (7) 
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We note that in addition to hexagon currents, new variables, a = ±1, which reside in 
triangles, appear. These variables count the multiple ground states. Now we return to 
calculation of the partition function (Q. We substitute the asymptotic formula for the triple 
products of Bessel functions ([7]) into Eq. (|6]) and rewrite the product of exponentials over the 
triangles as an exponential of the sum over the triangles. We also use the Poisson summation 
formula. Finally, we arrive at the following expression for the partition function 



Z(/3) = £ /exp 

MR t )},{m(R h )} J 
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27ri£j(R h )Q(R h ) 



Rv, 



£ (J(R h )-J(R h + B h ))' 

" Rh,B h 



ndJ(R h ), (8) 



Rv, 



Q(R h ) = m(R h ) + W a(R h + A t ) (9) 

6 At 

Here A t runs over all six triangles surrounding each hexagon with the centers Rh, Bh are 
six vectors that connect the centers of nearest hexagons. Note that centers of hexagons form 
a triangular lattice which is dual to the hexagonal lattice. 

Now one can integrate the partition function (|8]) over the currents in hexagons, J(Rh). 
This results in the expression for the partition function of the 2D Coulomb gas with quasi- 
charges Q(Rh) (Eq- (§)) positioned on sites of the triangular lattice Rh- Quasi-charges are 
1/3-multiple, this corresponds to the 27r/3-multiplicity of vortex rotations. 

At zero temperature, the integration over J(Rh) in (J|) yields conservation conditions 
IlR h ^(Q(R-h))) i- e -> m an y ground state, the sum of chiralities of triangles surrounding each 
hexagon is a multiple of 3. The problem of counting ground states is mapped onto that of 
coloring of the hexagonal lattice |J which was solved exactly 0. The exact number of ground 
states, Zn, is equal to 1 .46099^' 3 , where N is the number of spins. A naive approximation 
that assumes that chiralities of triangles surrounding each hexagon are independent and 
equally probable gives a good estimate Z N « (11/8) N/3 = 1.375 N/3 for the number of 
the ground states. In this estimate we neglect correlations between chiralities of triangles 
surrounding neighboring hexagons and farther correlations. The effect of those correlations 
can be estimated as the inverse number of the nearest neighbors on the triangular lattice, 



1/6. 



At finite temperatures, excitations against ground states, vortices with nonzero quasi- 
charges Q, appear. Both integration over hexagon currents, J(Rh), and summation over 
chiralities of triangles surrounding each hexagon contribute to the probability of vortex 
formation. Therefore, in addition to the standard estimate of the probability of vortex 
formation in unfrustrated XY magnets, one has to find the contribution due to various 
chirality configurations. This accounts for the high degeneracy of ground states in the 
Kagome antiferromagnet. The lowest energy vortices have quasi-charges Q = ±1/3. States 
with the sum of chiralities of triangles surrounding a certain hexagon equal to ±1 and ±4 
contribute to formation of such Q = ±1/3 vortices. The number of such configurations, -Zi,v? 
differs from the number of ground states, Zpj by some numerical factor, w\. We estimate the 
factor W\ the same naive way as we estimated the number of ground states, i. e. we assume 
that chiralities ±1 have equal and independent probabilities, w\ ~ 21/22. The precision of 
this estimate is again of order 1/6. 

Thus, the low-temperature partition function is the sum over states with zero quasi- 
charge and over states with quasi-charges equal to ±1/3. It is convenient to redefine \l/(Rh) = 
J(Rh)/(3iT) and K = (3/12. In terms of new fields \I/(Rh), the partition function reads 



This expression coincides with that obtained by JKKN ||. The quantity W\ plays the role 
of magnetic field yo in Ref. ||. In our case, the lattice Rh is dual to the lattice of hexagons, 
which corresponds to the initial hexagonal lattice. Hence, the problem of the KT transition 
on the Kagome lattice with the antiferromagnetic interaction is mapped onto that of the 
KT transition in a ferromagnet on the hexagonal lattice. 

In order to find the transition temperature and to estimate the role of chiralities we 
extract the short-distance part of the Green's function in the hexagonal ferromagnet (P~0|) . 
In the continuous limit, we obtain 




x 1 + 2w 1 ^cos(27r J Fr#(R h )) JJd*(R h ). 



(10) 




S 



Z= /D*(r)e /d2r L 2 ^ y '~ — (n) 



where h = 2wie~ Kn2 ^ 2 = 2w 1 e _ ' 37r2//24 . At the KT transition temperature, this is a small 
field, therefore it can be neglected according to reasoning by JKKN. From the usual renor- 
malization we find that the KT transition on the hexagonal lattice occurs when 

\/3/(2iQ = tt/2, i. e. 

TJJS 2 = V3n/72 = 0.0756. (12) 

Note that recent Monte Carlo simulations of the KT transition in the Kagome antifer- 
romagnet @] yield T C /JS 2 ~ 0.078. This is in a very good agreement with our exact result 
(0)- 

The existence of a new set of variables, chiralities, qualitatively changes the spin correla- 
tion function compared to that in "normal" XY magnets. Returning to the initial formula- 
tion of the problem (0), we consider the correlation function /C(r ) = (exp (i6(0) — i9(r ))). 
In terms of the integer- valued variables, n(r, a) (see (|3])), we arrive at an expression that 
differs from (||]), only by arguments of the 5-functions. Namely, for sites and r we get 

5>(0,a) = -5>(r ,a) = l. (13) 

a a 

instead of the conservation condition (^). This condition is equivalent to a pattern of currents 
which is a superposition of currents J(Rh), that flow in the Kagome lattice and obey the 
condition (§), and an unit current which is created in the point and is annihilated in the 
point r . Thus, the correlation function /C(r ) has the form 

Hr ) = ^e n 

^V^V {J(R)} (R^R*,A^A*) 

-fj(R+A)+J(R)(-/3) II Ij(R*+A*)+J(R*)+l(-/3), (14) 
(R*,A») 

analogous to that of Eq. (^). Here (R*, A*) are sites and vectors of the dual lattice such that 
A* crosses the path (0, r ) on the initial kagome lattice. Now we integrate over currents in 
the triangles in Eq. (|14|) as we did it before using the asymptotic formula (0). Instead of 
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Eq. (D we get /C(r ) = Z{(3, r )/Z(J3) where Z(f3) is given by Eq. (§) and Z(/3,r ) differs 
from by the additional contribution from the unit current running along the path 

(0,r ). 

The vertex contribution in the large tq asymptotics of the spin correlation function 
/C(r ) below the KT-transition point is negligible because the renormalization-group flow 
at T < T c yields that the effective constant h in (|ll|) is equal to zero. This is equivalent 
to the neglect of the first term in square brackets in Eq. (§) in the expression for Z((3,tq). 
Neglecting constraints on chiralities of triangles as we did before we immediately get a factor 
(cos27r/3) r °/ a = (— l) r o/ a 2 -r °/ a in the correlation function, where a is the Kagome lattice 
constant. The integration over J(Rh) in the r — > oo limit can be done in the spin-wave 
approximation and the result coincides with that by JKKN ||. Thus, in the low-temperature 
phase T < T c in the long-distance limit r /a^> 1 the spin correlation function reads 

/C(r ) oc (-l) ro/a 2- ro/a (r /a)- r/36Tc . (15) 

It decays exponentially with distance. The true order parameter of the KT-transition is 
the cubed spin ||, r/(r) = exp(3i0(r)). The correlation function of this order parameter at 
T < T c decays as a power of distance 

< 77(0)77(1-) >~ (r /a)" T/4Tc . (16) 

In conclusion, it is shown that the XY antiferromagnet on the two-dimensional Kagome 
lattice exhibits a Kosterlitz-Thouless transition and the transition temperature, T c , is eval- 
uated (Eq. (|12"D). It is found that the spin correlation function decays exponentially with 
distance even below the transition temperature (Eq. ( |I5|) ). Nevertheless, the order parame- 
ter correlations decay as a power of the distance in the low-temperature phase (Eq. (|T6"D). 
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